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Abstract 



q 

Analysis of electromagnetic planewave propagation in a medium which is a spatiotemporally homoge- 
neous, temporally nonlocal, isotropic, chiral medium in a co-moving frame of reference shows that the 
medium is both spatially and temporally nonlocal with respect to all non-co-moving inertial frames of ref- 
erence. Using the Lorentz transformations of electric and magnetic fields, we show that plane waves which 
have positive phase velocity in the co-moving frame of reference can have negative phase velocity in cer- 
tain non-co-moving frames of reference. Similarly, plane waves which have negative phase velocity in the 
co-moving frame can have positive phase velocity in certain non-co-moving frames. 
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1 Introduction 



Analysis of planewave propagation in a frame of reference that is uniformly moving with respect to the 
medium of propagation can lead to the emergence and understanding of new phenomcnons. In this respect, 
the Minkowski constitutive relations, as widely described in standard books [1, 2], are strictly appropriate 
to instantaneously responding mediums only [3]. For realistic material mediums, recourse should be taken 
to the Lorentz transformation of electromagnetic field phasors [4, 5]. 

In an earlier study on plane waves in a medium that is spatiotemporally homogeneous, temporally 
nonlocal, isotropic and chiral in a co-moving frame of reference, we reported that planewave propagation 
with negative phase velocity (NPV) is possible with respect to a non-co-moving frame of reference, even 
though the medium does not support NPV propagation in the co-moving frame [6]. That study applies 
strictly only at low translational speeds. In this paper, we demonstrate by means of an analysis based on 
the Lorentz-transformed electromagnetic fields in the non-co-moving frame, that our conclusion remains 
qualitatively valid for realistic mediums even at high translational speeds. 
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2 Planewave analysis 



We consider a spatiotemporally homogeneous, spatially local, temporally nonlocal, isotropic chiral medium, 
characterized in the frequency domain by the Tellegen constitutive relations [7] 



D' = e e' r E' + i^/i^'H' 



B' 



(1) 



in an inertial frame of reference £'. The relative permittivity e' r , relative permeability fj/ r and chirality 
parameter £' are complex-valued functions of the angular frequency u>' if the medium is dissipative, and 
real-valued if it is nondissipative [1, p. 71]; e and fi are the permittivity and permeability of free space, 
respectively. The electromagnetic field phasors are related by the Maxwell curl postulates as 



V x H' + iu'Tt' 

V x E' - iu'B' -- 

Our attention is focused on a plane wave, described by the £' phasors 

E' = E(, exp [i (k' • r' - u't')] 1 
H' = H cxp [i (k' • r' - Lu't')} J ' 



(2) 



(3) 



which propagates in the medium characterized by (1), with wavevector k' = kl k' and wavenumber k' = k k' r . 
There are four possibilities for the relative wavenumber: k' r € {k' rl , k' r2 , k' r3 , k' r4 } where 

Kl — y / e 'r^'r + C 

k' r2 — yje' r ^' r — C 



K r3 

y 



-V 



(4) 



With respect to frame £', the plane wave is assumed to be uniform; i.e., k' £ K 3 . 

Suppose that the inertial frame £' is moving at constant velocity v = wv relative to another inertial 
frame £. The electromagnetic field phasors in £ are related to those in £' by the Lorentz transformations 

[4,5] 

E = (E' • v) v + 7 [(I - vv) • E' - v x B'] 

"/ -\ v x E'" 

(1-vv) 'B> + — — 



B = (B' • v) v + 7 

H = (H' • v) v + 7 [(I - vv) • H' + v x D'] 

v x H' 



D = (D' • v) v + 7 



(Z-vv) »D' 



(5) 



where 7 = xx + yy + zz is the 3x3 identity dyadic, 7 = 1/ \/l — /3 2 and the relative translational speed 
(3 = v/c , with c = l/^/eoMo being the speed of light in free space. In terms of the £ phasors, the plane 
wave is described by 

E = E exp [i (k • r - uit)] 
H = H exp [i (k • r — ut)] 

The phasor amplitude vectors {Eo,Ho} and {Eq,Hq} are related via the transformations (5), whereas [4] 

r= [Z+(7-l)vv] T' + 7vt', (7) 
v • r 



(6) 



t = 7 K + 



k = 7 ^k' • v + — ) v + ( I - v v) • k' 
w = 7 (w' + k' • v) . 



(8) 

(9) 
(10) 
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Since k' e C for a dissipative medium, we have from (9) that k = k R k R + ikikj with fc/j e M, fc/ € R, 
k/? e R 3 , and k/ € R 3 , but k/j ^ kj in general; i.e., the plane wave is generally nonuniform with respect to 
E. Similarly, from (10) we have that u = uo R + iuoi with uo R <E R and ujj e M. Expressing the E phasors as 



E = {E exp [- (kj • r - w/i)]} exp [i (k/f • r - uigt)] 
H = {H exp [- (k/ • r - w/i)]} exp [i (k R • r - 



(11) 



we note that that the periodic propagation of phase is governed by k R and u>r, whereas attenuation or growth 
of the wave amplitude is governed by k/ and wj. By writing the phasor amplitudes as E = E ,r + «E / and 
H = H _r + iHoi, the corresponding cycle-averaged Poynting vector may be expressed as 



exp (-2k 7 • r) — x 

ZTT 



Eq/j x Ho/i 



rto + i 
J t 



cos 2 (k_R • r — ujRt) exp (2w/t) dt 



(E / x H 0R + E r x H /) x 

r - Lu R t) sin (k fl • r 



J t 



cos (k fl 

to + f 



ujRt) exp (2w/t) 



+Eq/ x Ho/ 



sin 2 (k/j • r — u> R t) exp (2w/t) rft 



for a cycle beginning at time t = to- The phase velocity is given by 

v p = — kR. 



(12) 



(13) 



Whether the plane wave has positive phase velocity (PPV) or negative phase velocity (NPV) in the refer- 
ence frame E is determined by the sign of v p • P: positive for PPV and negative for NPV. Criterions for 
determining whether the phase velocity is positive or negative with respect to the reference frame E' are 
presented elsewhere [8, 9]. 



3 Numerical results and discussion 



For the sake of illustration, let us consider the cycle-averaged Poynting vector evaluated at the point r = 
with the temporal averaging starting from to = 0; i.e., 



lr=0 



exp 



/ 4ttu>/ \ _ i 



(I 



- ^/) E /? x H /j 



87ro;/|a;| 2 

—lorlui (Eo/ x Hq/j + Eq/j x Ho/) + u> R Eq/ x Ho/ 



(14) 



Without loss of generality, let us assume that the plane wave propagates along the z' Cartesian axis; i.e., 
k' = z'. It follows then from the Maxwell curl postulates (2) that E lies in the x'y' plane with 



y'.E -zx'-E for k' r = k' rl ,k' r4 

V V 

"V2i "V3 



y • E = -iit' • E' for k 



Further, we take the velocity v to lie in the x'z' Cartesian plane as per 

v = x' sin 9 + z' cos 9 . 



(15) 



(16) 



In Figure 1 the distributions of PPV and NPV in the reference frame E for the dissipative scenario 
wherein e' r = 6.5 + 1.5i, £' = 1 + 0.2i and fi' r = 3 + 0.5i are displayed for k' r £ {k' rl ,k' r2 , k' r3 , k' r4 }. Clearly, 
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for all four values of k' r , propagation is of the PPV type for (3 = (i.e., with respect to the £' frame). As 
the relative translational speed increases, the phase velocity of the plane waves corresponding to k' rl and k' r2 
eventually becomes negative provided that tt/2 < 8 < tt. In contrast the plane waves corresponding to k' r3 
and k' r4 have NPV at sufficiently large values of (3 provided that < 8 < tt/2. 

Now, let us look at the scenario where the isotropic chiral medium supports NPV propagation in the 
co-moving reference frame. This situation arises when |Re {£'} | > Re {y^y^}, for example [8, 9]. We 
take e' r = 6.5 + 1.5i, £' = 10 + 2i and \J r = 3 + 0.5i. The corresponding distributions of NPV and PPV are 
mapped against [3 and 8 in Figure 2. We see that the plane waves corresponding to relative wavenumbers 
k' rl and k' r3 have NPV when (3 is small but have PPV when (3 is sufficiently large. On the other hand, the 
plane waves corresponding to relative wavenumbers k' r2 and k' r4 have PPV when (3 is small but have NPV 
when (3 is sufficiently large. 

To conclude, by means of the Lorentz-transformed electromagnetic fields, we have demonstrated that a 
plane wave with PPV in an isotropic chiral medium can have NPV when observed from a no-co-moving 
inertial reference frame. Similarly a NPV plane wave in the co-moving frame can be PPV from a non-co- 
moving frame. 
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Figure 1: The distribution of positive phase velocity (PPV) and negative phase velocity (NPV) in the E 
reference frame, in relation to [3 G [0, 1) and 9 e [0°, 180°), for k' r S {k' rl ,k' r2 , k' r3 , k' r4 }. Here, e' r = 6.5 + 1.5i, 
£' = 1 + 0.2i and n' r = 3 + 0.5i. 
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Figure 2: As Figure 1 but with £' = 10 + 2i. 
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